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ABSTRACT 

We compute the one loop self-mass- squared of a massless, minimally cou- 
pled scalar which is Yukawa-coupled to a massless Dirac fermion in a general 
conformally flat background. Dimensional regularization is employed and a 
fully renormalized result is obtained. For the special case of a locally de 
Sitter background our result is manifestly de Sitter invariant. By solving 
the effective field equations we show that the scalar mode functions acquire 
no significant one loop corrections. In particular, the phenomenon of super- 
adiabatic amplification is not affected. One consequence is that the scalar- 
catalyzed production of fermions during inflation should not be reduced by 
changes in the scalar sector before it has time to go to completion. 
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1 Introduction 



Essentially all quantum field theoretic effects can be understood through 
classical interactions of the virtual particles whose existence is required by 
the uncertainty principle. In general one expects quantum field theoretic 
effects to become stronger the longer virtual particles live and the more 
probable it is for them to emerge from the vacuum. For example, vacuum 
polarization arises due to the polarization of charged virtual particles in an 
external electric field. The largest effect derives from electron-positron pairs 
because they are the lightest charged particles and therefore live the longest. 
One can even understand the running of the electromagnetic force from the 
incomplete polarization of the longest wavelength virtual pairs in the field 
provided by two very close sources. 

The expansion of spacetime affects quantum field theory by lengthening 
the time virtual particles can exist, and sometimes by altering the probability 
with which they emerge from the vacuum. The first effect can be understood 
from the energy-time uncertainty principle. In co-moving coordinates the 
geometry of a homogeneous, isotropic and spatially flat universe is, 



Although t measures physical time, physical distance is Ax multiplied by 
the scale factor a(t). Because spatial translation invariance is still a good 

— * 

symmetry one can label particles by their wave vectors k as in flat space. 
However, the physical momentum of such a particle is k/a(t), and one can 
think of the corresponding energy as, 



The spontaneous emergence of a pair with wave numbers ±k at time t will 
not lead to a detectable violation of energy conservation provided the pair 
persists no longer than a time At defined by the equation, 



Hence we conclude that the expansion of spacetime always increases the time 
which a virtual pair can persist. 



ds 2 = —dt 2 + a 2 (t)dx ■ dx . 



(1) 




(2) 




(3) 
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Just as in flat space the persistence time At is longest for massless par- 
ticles. In an expanding geometry this means m <C H(t), where the Hubble 
and deceleration parameters are, 

*<H ■ f = -i-|- w 

In this case the equation for persistence time takes the form, 

..=.„ /"*+ At dt' 

m = =^ 2\\k\\ / -r-r = 1 . 5 

11 11 Jt a(t') v ; 

The integral in © is just the conformal time interval At], defined by the 
change of variables drj = dt/a(t). The dependence of Arj upon At is con- 
trolled by sign of q(t). For q(t) > (decelerating expansion) Ar] grows 
without bound; while for q(t) < (inflation) it approaches a finite constant. 
Hence we conclude that any sufficiently long wavelength virtual particles 
which are produced during inflation can persist forever! 

Whether or not effectively massless particles actually engender stronger 
quantum effects during inflation depends upon the probability with which 
they emerge from the vacuum. Almost all massless particles possess a sym- 
metry known as conformal invariance which means that physical processes 
are the same as in flat space when expressed in conformal coordinates, 

dt 

dt] = — — ==>- ds 2 = a 2 (^—dr] 2 + dx ■ dxj . (6) 

Hence the number of virtual particles which emerge from the vacuum per 
conformal time is the same as the constant flat space rate we might call T. 
It follows that the rate per physical time falls off, 

dN _drj dN T 
dt dt di] a(t) 

Hence we conclude that while any sufficiently long wavelength, massless vir- 
tual particles which happen to emerge from the vacuum can persist forever 
during inflation, very few conformally invariant particles will emerge. 

The two massless particles which are not conformally invariant are min- 
imally coupled scalars and gravitons. For these particles it turns out that 
the rate of emergence per unit physical time is unsuppressed, so we expect 
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quantum field theoretic effects from them to be stronger than in flat space. 
In fact it is quantum fluctuations of these fields which are responsible for the 
scalar pQ and tensor [2] perturbations predicted by inflation [3J H] . These are 
tree order effects. At two loop order it can be shown that a massless, mini- 
mally coupled scalar with a quartic self-interaction experiences violations of 
the weak energy condition on cosmological scales E] and that quantum 
gravitational back-reaction slows inflation [7j. 

Conformal invariance in the free theory need not rule out significant quan- 
tum corrections if the conformally invariant particle couples to one which is 
not conformally invariant. Two examples of this have been studied recently. 
In the first, electrodynamics — which is conformally invariant in D = 4 
spacetime dimensions — is coupled to a charged, massless, minimally cou- 
pled scalar. The one loop vacuum polarization [HJ E] induced by the latter 
causes super-horizon photons to behave, in some ways, as though they had 
nonzero mass JUj- This engenders no physical photon creation during in- 
flation but leads instead to a vast enhancement of the 0-point energy of 
super-horizon photons which may seed cosmic magnetic fields after the end 
of inflation [THIP2]. 



+ 



x 
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Fig. 1: One loop contributions to jEj (x;x'). 

The other example consists of a massless fermion which is Yukawa-coupled 
to a massless, minimally coupled scalar ^3]. Massless fermions are confor- 
mally invariant (in any dimension) so they are not much produced by them- 
selves. However, the fermion self-energy (Fig. 1), and the effect it has on 
the quantum-corrected fermion field equations, show that the massless, min- 
imally coupled scalar catalyzes the emission of scalar-fermion-anti-fermion 
from the vacuum. In this paper we study what effect the fermion has on the 
scalar through the one loop self-mass-squared (Fig. 2). This can hardly alter 
the non-conformal coupling but it might induce a nonzero scalar mass. If 
this mass became large enough sufficiently quickly it could cut off the scalar- 
induced fermion creation. We will show that this does not happen. The one 
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loop scalar self-mass-squared can be renormalized so that there is no signifi- 
cant change to the scalar mode functions, and higher loop corrections cannot 
become significant soon enough. 

In the next section we use the Yukawa Lagrangian to derive the relevant 
Feynman rules for an arbitrary scale factor a(t). In section 3 we compute 
the renormalized scalar self-mass- squared at one loop order. Although our 
result is valid for any a(t) we show that it reduces to a manifestly de Sitter 
invariant form for the locally de Sitter case of a(t) = e Ht . In section 4 we 
use the self-mass- squared to study one loop corrections to the scalar mode 
functions. The consequences of this result are discussed in section 5. 




Fig. 2: One loop contributions to M 2 (x;x'). 



2 Feynman Rules 

We begin by reviewing the conventions appropriate to Dirac fields in a non- 
trivial geometry. In order to facilitate dimensional regularization we make 
no assumption about the spacetime dimension D. The gamma matrices 7^ 
{b = 0, 1, . . . , D — 1) anti-commute in the usual way, {■y b , r y c } = —2rj bc I. 
One interpolates between local Lorentz indices (b, c,d, . . .) and vector in- 
dices (lower case Greek letters) with the vierbein field, e^{x). The met- 
ric is obtained by contracting two vierbeins with the Minkowski metric, 
9nA x ) = e fj.b(x)e uc (x)r] bc . The vierbein's vector index is raised and low- 
ered by the metric (e M 6 = g^ v e v ^) while the local Lorentz index is raised and 
lowered with the Minkowski metric (e b = r) e^). The spin connection and 
the Lorentz representation matrices are, 

A, bc = e\(e uc ,, - T Vpc) , J bc = 71 • (8) 
Let 4>(x) represent a real scalar field and let ipi(x) stand for a Dirac field. 
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In a general background metric the Lagrangian we wish to study would be, 

-\^4> 2 RV=g - h<frH>V=g , (9) 

where ip = ^7° is the usual Dirac adjoint, £0 is the bare conformal cou- 
pling and fo is the bare Yukawa coupling constant. Note that we do not 
require mass counterterms because mass is multiplicatively renormalized in 
dimensional regularization. 

The geometry of interest is the very special form associated with the 
homogeneous and isotropic element. By defining a new time coordinate dr] = 
dt/a(t), the metric of this geometry can be made conformal to the Minkowski 
metric, 

ds 2 = a 2 (^—dr] 2 + dx ■ dxj . (10) 

A convenient choice for the associated vierbein is, = ai]^. With these 
simplifications the spin connection assumes the form, 

e^b = ar]^ b =>• A^ cd = (r]^ c d d - r]^ d d^j ln(a) . (11) 

And our Lagrangian reduces to, 

+ ^ (D-l)(2a flu a D - 3 + (D - 4)a ^a^rT '4> 2 - / o a D 0^. (12) 

Renormalization is facilitated by introducing the renormalized fields, 

= ^/~Z<t>R an d ^ = \f~Z2ipR ■ (13) 
This brings the Lagrangian to the form, 

+ \zUD-l)a D ~ 2 (2 d — + {D-^—y R - VzZ 2 f a D <f>R^ R . (14) 

Note the Dirac slash notation (^ = 7^<9 M ) and the convention — used hence- 
forth — that indices are raised and lowered with the Lorentz metric (e.g., 
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qh = fj^dy). Note also that d = <9 M <9^. We now define the counterterms, 

Z = 1 + 5Z , Z 2 = 1 + 8Z 2 , (15) 
VzZ 2 f Q = f + 6f , Z£„ = + <J£. (16) 

Note that the conformal coupling enters only as a counterterm because we 
want the scalar to be minimally coupled. We can now express the Lagrangian 
in terms of primitive interactions and counterterms, 



d,.a d^a 



-Uza^d^nO^R + U^D~l)a D - 2 (2 — + (D-A) 
2 2 v a a a 

+^ 2 (aV^)^( a V^) _ SfaP^^R. (17) 

We shall need the counterterms on the second line; we will not need those 
on the third line. 

The diagrams of Fig. 2 do not require the scalar propagator. The fermion 
propagator can be determined by noting from the second term of (117(1 that 
the combination a~^ 1 ip R behaves like a massless Dirac field in flat space. It 
follows that the Feynman propagator of ip is just a conformal rescaling of the 
flat space result, 



iSj 



[x-x') = (aa')^S^ | r( £ # 1} [Aa^foa/)] 1 "' j , (18) 



2ttt [Ax 2 {x;x')} 



aa')^ - l ' * ■ (19) 



Here Ax^ = rj^ ll ,(x u — x' u ) and the conformal coordinate interval is, 

Ax 2 (x; x') = \\x-x 1 1| 2 - {\r)-r)'\ - iS) 2 . (20) 

Note that we label the spacetime position with the D-vector x^ = (r),x). 
The split index notation in i[iSj](x; x') indicates that the first index (i) trans- 
forms according to the local Lorentz group at the first coordinate argument 
(x^) whereas the second index (j) transforms at the second argument (x ,fl ). 
The interaction vertex derives from the —fa D 4> R ip R tp R term of (fTTj) . 

-ifa D 5 ij . (21) 
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We also require the scalar field strength renormalization and conformal coun- 
terterms, 

i5ZdJa D - 2 d fl 8 D (x - x')) , (22) 



and, 

l 5i(D-l)a D - 2 (2— + (D- 4)^— )d D (x - x') . (23) 



a a a 

Let us note, for future reference, that the choice 



makes the two counterterms sum to a simple form, 
iSZdJa^S^x -x')) + i5Za D - 2 (^--l) j — 

2 ' a a 



+ 2)^—}5 u (x-x')=i5Z(aa')-- 1 d 2 5 v (x-x') . (25) 

V 9 / n n 



Because the "in" (t — > — oo) vacuum is not equal to the "out" (t— » +oo) 
vacuum in this background it is desirable to compute true expectation values 
rather than in-out matrix elements. This can be done covariantly using 
a simple extension of the Feynman rules known as the Schwinger-Keldysh 
formalism ^3J El QE1 E! • Briefly, the end of each line has a polarity which 
can be "+" or "— ". Vertices are either all + or all — . A + vertex is the 
familiar one from the Feynman rules whereas the — vertex is its negative. 

Propagators can be ++, H — , — h or . Each propagator can be obtained 

from the Feynman propagator by replacing the conformal coordinate interval, 
Ax (x\ x'), with the interval of appropriate polarization, 

Ax 2 ++ (x; x') = \\x — x\\ 2 — (\rj — r]'\ — i5) 2 , (26) 

Ax 2 + _(x;x') = \\x-x f - (r]-r]' + i5) 2 , (27) 

Ax 2 _ + (x;x') = \\x-x'\\ 2 - {ji-i -i5) 2 , (28) 

Ax 2 Ji;x') = \\x-x f - (\r)-r]'\ + i5) 2 . (29) 

External lines can be either + or — in the Schwinger-Keldysh formalism. 
Hence every N-point 1-particle-irreducible (1PI) function of the in-out for- 
malism gives rise to 2 N 1PI functions in the Schwinger-Keldysh formalism. 
For every field <p(x) of an in-out effective action, a Schwinger-Keldysh effec- 
tive action must depend upon two fields — call them <fi+(x) and <f)-(x) — in 
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order to access the appropriate 1PI function [T%1 ITT^ I2*U| . If external fermions 
are suppressed, the effective action for our model takes the form, 



0_] = S[<f> + ] - S[<j>.) 



X - Jd 4 x Jd 4 x' 



x 



lx' 



+ 0(cf>l), (30) 



<p + (x)Ml + (x; x')(p + (x') + <fr + (x)MlXx; x')(, 

(x)M 2 + (x; x')(j) + {x') + 0_(x)M 2 _(x; x')<p_(x' 

where S is the free scalar action. The effective field equations are obtained 
by varying with respect to <p + and then setting both fields equal (HI UH 120] , 



\x 



<t>±-- 



^(a 2 ^0(; 



x))- <Fx' M; + (x;x') + M;_(x;x') 0(x') + O(0 2 



(31) 

It follows that the two 1PI 2-point functions we need are M 2 + (x; x') and 
M 2 _(x; x'). Their sum in ()31|) gives effective field equations which are causal 
in the sense that the two 1PI functions cancel unless x ,fl lies on or within the 
past light-cone of x^. Their sum is also real, which neither 1PI function is 
separately. 



3 Renormalized One Loop Self-Mass 

In this section we compute and fully renormalize the scalar self-mass-squared 
at one loop order. Our result applies for any scale factor a(t). For the special 
case of de Sitter (a(t) = e Ht , with constant H) we give a manifestly de Sitter 
invariant form for M^ + (x; x') and M^_(x; x'). 

Using the Feynman rules of the previous section we see that the ++ and 
H — contributions from the first diagram in Fig. 2 are, 



^fa D 



i,Sj 



£c; x')(rif( 



ID 



■i Si 



[x ; x) 



±+ 



/ 2 r 2 (f) 



(HI 



IT 



D Ax 2 ^ 



(32) 

where a is the scale factor at conformal time r\ and a! at rj'. One cannot 
yet take the spacetime dimension to D = 4 because this term is too singular 
to give a well-defined integral in (|31)1. To make it less singular we extract 
derivatives with respect to the un-integrated coordinate x M , which can be 
pulled outside the integral. The key identity can be stated without regard 
to ± variations, 

1 1 d 2 (^-^). (33) 



Ax 



2a 



4(a-l)(a-f) 



Ax 



2(a-l) 
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Two applications of this identity give us, 

T mt) aa' /Wr'(f-l) u 1 x 

+ n D Ax 2 + ^-V + lQir D (D-3)(D-4) ^Axl^- 3 " ' 1 f 

At this point we could take the limit D = 4 were it not for the explicit factor 
of (D—4) in the denominator. 

It is now necessary to distinguish the ++ case, which has a one loop 
counterterm, and H — case, which is free of primitive divergences. The trick 
for obtaining the renormalized result in each case involves adding zero using 
the identities, 

1 \ «47TT , s , „o/ 1 



nA^) = rtf^T) 5 {x ~ x>) and 9 "a^) =0 - (35) 

We can therefore write the two self-mass-squared's as, 

■ M 2 ( /w r 2 (f-i) f i 

—lM*\x\ X = — ; - = — -f; — —O 



16tt° (Z>-3)(D-4) ^A4 ( + D - 3) 



+i 6^D-l)a D - 2 (2— + (D- 4)^—)5 D (x - x') + 0(f), (36) 



Note the appearance of the dimensional regularization mass scale \x. 

By comparing the primitive divergence in (}36|) with the simple countert- 
erm (f23|) that results from the relation (|24j) . we settle on the following choice 
of counterterms, 

^ = -A^rMz + + 0(f) and SZ = f -^l r /| Z 2 ) + (/ 4 ). 



(38) 

Here <5£ fln is a finite, order / 2 contribution we shall fix later. With this choice 
the ++ result becomes, 
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16tt d (£>-3)(£>-4) \Axl { r- V Ax? + T ~ i; 

if 2 ^(~2~^) D-4 I I I a2xD/ A 

r (A o\ I 2 [aa-{aa)2 )d 5 (x-x) 

+i5^ n 6ad 2 a5\x-x') +0{f) . (39) 

Even though the bare scalar is minimally coupled, the divergent parts of the 
one loop counterterms are interpretable as the field strength renormalization 
of the conformal kinetic operator! We might have anticipated this from the 
fact that only fermion propagators enter the primitive diagram at one loop 
order, and massless fermions are conformally invariant in any dimension. 
Higher loop diagrams such as those of Fig. 3 involve the scalar propagator, 
which breaks conformal invariance, so we do not expect the conformal relation 
(j23| between and 5Z to persist at higher loops. 

At this stage we take the limit D = 4 facilitated by the identities, 

D-4 , n . W,,2A~2> 



fi D ~ 4 (D xln(/i 2 Ax 



£-2)^™ + 0((Z>-4)>), (40) 



Ax 2(D-3) Ax 2(f-1) V2 J Ax 

aa'-iaa')^- 1 = -(—-2)aa! ln(aa') + 0((L>-4) 2 ) . (41) 

The factor of ln(aa') in the second relation is reminiscent of the non-local 
conformal anomaly [21] and derives from the same dimensional mismatch 
between primitive divergence and counterterm. Putting everything together 
and taking the limit D = 4 gives, 

Ml(x;x') = ^4d^ Hfi2Ax ^ 



' +y ' 1 32vr 4 V Ax 2 ++ J 

J ln{aa')d 2 S\x - x') - 5£ &l fiad 2 a5\x - x') + 0(/ 4 ) , (42) 



8tt 2 

Ki^) = -^(^^) + 0(f). (43) 

We now specialize to de Sitter background, i.e. airf) = —l/Hrj with H 
constant, to show that the self-mass- squared can be expressed in a manifestly 



10 



de Sitter invariant form. The de Sitter invariant, conformal d'Alembertian 
is. 



6 



-gR 



ad 2 a 



A simple function of the invariant length i(x\ x') is, 

y(x\x') = Asm 2 (^-H£(x; x')^j = aa'H 2 Ax 2 (x; x') 



(44) 



(45) 



The polarized forms of this length function are obtained by simply replacing 
the coordinate interval, Ax 2 (x;x'), with the interval of appropriate polariza- 
tion from (J26j) through (|29|) . Using these invariants, (jl2|) and (}43|) can be 
rewritten as 



M 2 ++ (x; x') 



Ml_(x; x') 



32tt 4 

j_PH 2 
32tt 



-6^Ry/=g-6\x - x') + 0(f), (46) 



4 W 



, \n[y + _(x;x')^ 2 /H 2 } 



y+-{x;x' 



+ 0(f 



(47) 




4 Effective Field Equations 

In this section we substitute our results (j42H43|) for M 2 ± (x; x') into the ef- 
fective field equation (|3*T]) and work out the result for a spatial plane wave. 
Most of the analysis is valid for arbitrary scale factor a(t). Only at the end 
do we specialize to the locally de Sitter case and make a choice for <5£ fln which 
keeps corrections to the wave functions small at one loop order. 
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The linearized, effective field equation is, 
f 2 a 

d^(a 2 d^<p(x)) + J — [\n(a)d 2 (a(f)(x)) + d 2 (\n(a)a(f)(x))] + 5^ n Qa(d 2 a)(f)(x 

if 2 a 
32tt 4< 



& i l<rx'0{rf-rn)a'<l>{x!) 



ln(// 2 Ax 2 ) ln(/x 2 Ax^ 



Ax 



++ 



Ax 2 



+ 0(/ 4 ) = 0. (48) 



Here rjj is the initial conformal time (corresponding to t — 0) at which the 
state is released in free Bunch-Davies vacuum. The first step in simplying 
this equation is to extract another d'Alembertian from the nonlocal term in 
square brackets, 



ln(/i 2 Ax 2 + ) ln(// 2 AxJ 



Ax 



Ax z 



d 2 r ln 2 ( M 2 Ax 2 _ + )-21n( At 2 Ax 2 + ) 
-ln 2 (/i 2 Ax 2 ) + 21n( y u 2 Ax 2 



We now define the coordinate intervals Arj = rj — rj' and Ax 
recall the ++ and H — intervals, 



• (49) 
x — x'll and 



Ax 2 + = Ax 2 - (\Ar]\ -iSf 



and 



Ax 2 _ = Ax 2 - (Arj + i8f . (50) 



When rf > rj we have Ax 2 + = Ax 2 _, so the ++ and H — terms in (|49|) 
cancel. When rf < rj and Ax > At] (spacelike separation) the arguments 
of the logarithms become positive, real numbers for 5 — > 0, so there is also 
cancellation. Only for rf < 77 and Ax < At] (timelike separation) do we 
acquire a nonzero result through the relation, 



9(Ar) - Ax) ln(/x 2 Ax 



9(Ar] - Ax) I In fi 2 (Arj 2 - Ax 



± in 



(51) 



Hence the square bracketed term in ()48|) can be written as, 
ln(/i 2 Ax 2 + ) ln(/i 2 Ax 2 _; 



Ax 2 + Ax 2 . 



y<9 2 <| 6(Ari-Ax)(\n[fi 2 (Ari 2 -A 



(52) 



Substituting this relation into (jjHj) gives the manifestly real and causal equa- 
tion, 

f 2 a 

9 M (a 2 ^0(x)) + J — [ln{a)d 2 (a<f){x)) + d 2 (ln{a)a<f){x))] + 5£ Rn 6a{d 2 a)<f){x) 



2 6 7T 3 



^fdrfairf) f rf 3 x'0(V,x')(ln[^ 2 (Ar/ 2 -Ax 2 )]-l) + 0(f) = 0. (53) 
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Because the background geometry is homogeneous, isotropic and spatially 
flat, we can build up an arbitrary solution from a superposition of spatial 
plane waves of the form, 

<f>(ri, x) = g(v, k)e ik ' s . (54) 

Evaluating the derivatives of this in the first two terms of (|53|) is straight- 
forward. The nonlocal term, involving the integral, is more complicated. To 
begin we make the change of variable y = x'—x and extract the spatial phase 
factor, 

Pa 



2 6 7T 3 



d 6 [dv'a(ri')g(ri',k) f d 3 x r e*' s ' (ln[// 2 ( A77 2 - Ax 2 )] - 1 

Jr/i J Ax<Ar) ^ 



P a Jk-xi 



2 6 7T 3 



e^(d 2 +k 2 ) 3 d V 'a(v')9(v',k) d 3 ye^ ln[/i 2 (A^ 2 -y 2 )]-l .(55) 



s/1I<Aj? 



We next perform the angular integrations and make the change of variables 

y = A772;, 

2 6 7T 3 



d 6 [ V dr ] , a(r] , )g(T] , ,k) [ rfVe^'(ln[/i 2 (Ar/ 2 -Ax 2 )]-l 

Jrji JAx<Arj ^ 



e ik -*(d 2 + k 2 ) 3 rdr]'a(ri')g(ri',k)l dy»- sin(%) ln[// 2 (Ar/ 2 -t/ 2 )], (56) 



fa 
2% 2 

P a Ak-ZfcO, 1 i„2\3 rV 



y . 



2% 2 



k 



x / dz At/ 2 - sm(kAr]z) [ln(/i 2 Ar/ 2 ) + ln(l - z 2 ) - 1 



e lk - x (d 2 + k 2 Y 1drfa{rf)g{rf, k) 



The integral over z is facilitated by the special function, 

£(a) = / dz zsm(az)\n(l — z 2 ) 
Jo 

si(2a) + - 



2 ■ 1 r ■ 1 

— sm a cos a + a sm a 



H — - sin a — a cos at] 



ci(2a) — 7 — In 



a 



(57) 



(58) 



(59) 



Here 7 is the Euler-Mascheroni constant and si(x) and ci(x) stand for the 
sine integral and cosine integral functions, 

sint 



si IX) 



Cl(X) 



00 , sint 

dt 

t 



71 



h / dt 



t 

r°° , cost , r x , 

- I dt = 7 + In x + / dt 



t 







cos t — 1 



(60) 
(61) 
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Making use of these relations and performing the elementary integrals gives, 

f 2(X «6 Hi, 



■^d 6 dri'a(v')g(v',k) / d 3 x'e kS> (ln^Ary 2 -Ax 2 )]-1 

2 7T Jj^ J Ax<Ar/ ' 



ik-x 

2% 2 



drj'a(r]')g(r]', k) 



x 



(sin(fcAr/)-A;Ar7cos(A;Ar7))(21n(/iAr/)-l) + (kAr]) 2 £(kAr)) 



(62) 



The next step is acting the derivatives. Because the integrand in (|p^j) 
vanishes at 77' = 77 like ln(Ar/)Ar/ 3 the first three derivatives commute with 
the upper limit, 

f 2(l 



4-^a 6 /Va(V)^(V,fc) / rfVe^'fln^^A^-Ax 2 )]-! 

2 b 7T d JAx<Ar) V 

r ' 1 e ifc ^ (<9 2 + /c 2 ) (flb + ik)(da - ifc) [dr)'a(r)')g(r}', k) 



2H 2 k 
x< sin(A;A?7) 



Ill 



'2/i 2 A^ 



+ ci(2A;A?])-7 



7T 



-cos(A;Ar7)(si(2A;Ar/) + - )} ,(63) 



A. 

2% 2 



e^ 5 (<9 2 + A; 2 ) (5b + zife) fdrfairf) g(rj', k) e~ ikAri 



2 ln^/xAr?) + jj^tfe ( ^— ^ ) | • (64) 



The term containing the logarithm in ()64|) is divergent for 7/ = rj. We isolate 
the divergence using, 

ln(2/iA?7) = ln(-2/ir/) + ln(l - ^) . (65) 
We can now act the operator (do+ik) on the non-singular terms in (JMj) . 

2 6 7T 3 



4tt 2 



d % \ r drja{rj)g{rj,k) f d 3 x' e^ 3 ' (\n\ji 2 {Arj 2 - Ax 2 )} -1 

Jr?i J Ax<Arj ^ 

a(v)g(v, k ) ln(-2/xr/) + z JJ dr]'a(ri')g(r}', k)± sin(A;A?7)| 
+ («% + zfc) JJ dr)'a(r)')g( V ', k)e~ ikA n In (l - J 



x(«9 2 + £; 2 



. (66) 



Had we instead acted (<9 +z/c) before (d — ik) the result would be, 
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( B 2 ,2) / a (v)9{v, k) ln(-2/i?7) - i J* dr]'a{<n')g(r}', k) ^ sin(fcAr?) 1 
X V° + ;\ +(d -ik)^d V 'a( V ')g( V ',k)e^nn(l-^) j' { j 

A simpler expression results from adding half of (|66|) with half of (J67|) . 

■£±&rdrfa{jf)g{j(,k) f d z x'e^ (\n[^{Arf-Ax 2 )\ -l) 

r a(f])g(i], k) ln(-2//77) ^ 

= -£i| e ** («9 2 + A; 2 ) +<% JJ d V 'a( V ')g(rj', k) cos(A;Ar ? ) In (l - ^) . ( 68 ) 
1 " " [ d V 'a( V ')g( V r , k) sm(kA V ) In (l - J 

The top line of (J55|) is comparable to the local one loop terms in ([53)1 . 
Extracting it is what we have worked so hard to do. The remaining, nonlocal 
terms make only small contributions at late times. They can be simplified 
by first expanding the trigonometric functions using the angular addition 
formulae, then partially integrating to shield the logarithmic singularity, and 
finally bringing another derivative inside, 

n f Art \ fw / An 

d / drj'a{rj )g{rj , k) cos(fcA^) ln( — -J + k / dr\ a\r\)g[r\ , k) sm(kAri) ln( — - 

Jr\i TJ J rji TJ 

rv / An \ 

= cos(krj)do drfa{r]')g{ri ,k) cos(kr]') ln( — -J 

Jrii v — TJ 

+ sin(kr])dof drf a(rf)g{rf , k) sm(krj') lnf — (69) 

Jlji V — TJ 

1 )—a(r)i) cos(k(n - r]i))g(r)i, k) 

Vv V v ' 

+~ jTVln(^) A {7/0(7/) cos^At?)^', &)}. (70) 

It is now time to combine terms and give the general result for the lin- 
earized effective field equation (JHHj) specialized to a spatial plane wave (|54jl . 
We denote conformal time derivatives with a prime, 

aJa^fcfaAOe**)! = - a *e iS *(gr + 2-g t + k*g), (71) 



a 
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-(d 2 + k 2 )(ag) = g" + 2-g' + k 2 g, 

CI (X 



(72) 



(d 2 + k 2 ) (\n(-2m)ag) = (d 2 + k 2 )(ag) + 2 -{ag)' - ±ag . (73) 
Combining everything, and deleting an overall factor of —a 2 e lk ' x , we obtain, 
= g" + 2-sf + k 2 g + ^ ln(-2^a) [g" + 2-g' + (- + k 2 )g 



An 2 



,a! 2x . ( a! 2 a" 2a! 1 

(- + -)9 +{7^ + 7r + 2 



/ 2 



47T 2 ( 



a 



• (><*£=..—// 
la rja 7f J J a 

ln(l - f a(r/i) cos(k(r]-r]ifjg(r] u k) 

+ H drf ln(^) £ [^a(r/') coB(kAr,)g(rf, fc) 



(74) 
+ 0(/ 4 ). 



One can better infer asymptotic behaviors if an extra factor of a 2 is ex- 
tracted and the equation is converted to physical time t, 



k 2 

g + 3Hg+-g + ^ ln(-2^a) [5 + ZHg + (2# 2 + if + ^)g 



f 



./ 



4tt 2 

2 \ . /3 tt9 3 • 2iJ 
# + — )g+(-H 2 + -H+ — 
r/a / v 2 4 77a 



^(^ +3 4 +2 ^ +ff 4){ in ( i -^^ cos (*("-"-»^ fc » 



r*' in (^)^[$|cos(^) 9W ,*)]} +0 (/')=o. 



(75) 



It is also useful to recall the slow roll expansion for the conformal time, 

alt' -1 



V 



-jT 



a(t') H(t)a(t) 



(76) 



Hence the combination r/a(t) is only slowly varying during inflation. The 
tree order mode function approaches a constant at late times so we will 
keep one loop corrections small by making finite renormalizations to cancel 
any undifferentiated mode functions which are not suppressed by slow roll 
parameters. The best choice for this is, 



fi = l -H{U) = ±Hi 



and ££ fln 



P 



32tt 2 



(77) 
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With these choices the effective mode equation becomes, 



k 2 



g + 3Hg+^g + £j In(-^a) [g + 3Hg + (2H 2 + H + ^)g 



4vr 2 



H+^)g+(3H 2 + 3 -H + ^-4- 2 
v 4 rja r] z a z 



rja 



f 2 (& +?>H ® + 2H 2 +H+K 



4tt 2 \dt 2 dt 



I ln(l- — ) — cos( k{r]-i] i ))g(r] i , k) 
I v rji' rja v ' 



Ar] ^ d \r}'a{t') 
i]a(t) 



cos(kAr])g(r)', k) + 0(f 



0. 



(78) 



The next step is to solve the equation perturbatively, 

P 

g{v,k) = g (v, k ) + Y^9i(v, k ) + ■■■ 

Because the tree order mode function obeys, 



(79) 



k 2 

g + 3Hg + — g = 



(80) 



the equation for the one loop correction is, 



g 1 + 3Hg 1 + ^g 1 = -\n(-^r ] a) (2H 2 +#) # - (#+ A) «? 

1 



(3H 2 + 3 -H^ 2/7 



r]a r] 2 a 2 



go 



d 



k 2 



— + 3H— + 2H Z + H+ — 



x 



dt 2 dt 
\n(l-%)^cos(k( V - Vi ))g ( Vi ,k) 

cos(kAr))g (r]',k) 



V'a(t') 
r)a(t) 



(81) 



Although the solution to (JHU|) has been obtained for a general scale factor a(t) 
[22] the expression is too complicated for the integral in (jHTjl to be evaluated 
in closed form. For the special case of de Sitter background (a(t) = e Ht = 
— 1/Hi], with H constant) the tree order mode function is, 



H 

7W 



ik 
~Ra~ 



exp 



ik 
~Ha~ 



H 

7W 



(l + ikr] 



-ikr] 



(82) 
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Of course many other terms in the effective mode equation (J75j) also simplify 
in de Sitter background, 



k 2 

g\ + 3Hg^ — -g x = Hg 
a 1 

( d 2 d 

■ + 3H— + 2H 2 + - 



hi(l- cos(k(r]-T] i ))go(Vi, k) 



dt 2 ' ""dt ' — ' a 2 J \+ j* dt> ln(^) § [cos{kAr,)g {i , k) 
To begin evaluating the nonlocal term we note the differential identity 



83) 



ln (z^) jnfa&^vMrf, k)} = I In(^) cos(kArj)g (T]', k) 







Arj- 



drj' 



+ 



V 

I 1 

T)' Af) 



cos(kArj)g (v',k) . (84) 



The 1/rj' term on the right hand side of (|84|) can be expressed as, 



1 H 

— cos( k An) g ( ri' , k) = ■ 





-i2krf _i 


rf 





+ 



d_ 

drf 



2 cos(kri) In(-Hrf) - - e ik ^~ 2 ^ + ikrfe~ ikr) 



(85) 



A similar expression can be obtained for the 1/Arj term on the right hand 
side of 



1 H 

—— cos( k An) gdrf, k) = 

A77 v /;yuv " ; (2k) 



[l+iki])e 



-ikrj 



e i2kArj_^ 



A77 



+ 



d_ 

drf 



-2(l + ikri)e- ikr >ln(HAri) + -e ik ^~ 2 V) - ikrj 'e~ ikr <] . (86) 



It follows that the bracketed term in (1831 is 



ln(l -J cos(k(i] -11^ g^iii, k) + J drf\n(— ) — [cos(kAr])g (rf ,k) 



rf ' drf 

ln(l + Hr])g (r], k) + \n(-Hrj) cos(ki])g (0, k) - g (ij, k) 



--g o (0,k)e ikr > / dz 



-2kr]i 



2kt] 



e lz -l 



2kAr) i 

+ ^9o(v,k) I dz 



e lz - 1 



• (87) 
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The integrals in this expression could be written in terms of the sine and 
cosine integrals (|6UH61|) but it is simpler not to do this. 

It remains to substitute the nonlocal source term (|57j) into the equation 
()83|) for the one loop mode function. Acting the derivatives gives a compli- 
cated expression which considerable effort brings to the form, 

k 2 n ( ( p 2ikAr !iA_l \ ~] 

H~% + 3tf-V + = (g* - 9o ) Ha) + H-% jl-f - _ f j \ 



+g { -2ikr] - 2 ln(l + Hrj) 



'l-ikri)e 2ikAr >> + l (e 2ikAr >* + l 



+ 2 9 ° 



-2k Vi re iZ -l 
dz 

2kr) 



a — 1 
1 



+ 



2(a - V 



:9o 



2kA Vi re iZ - 1 
dz 



All the mode functions in this expression are evaluated at rj. From the 
asymptotic late time expansion for the mode function, 



9o(v, k ) 



H 



'2k z 



1 / k \2 i / k \3 



(a*) 



(89) 



we see that the leading form of the late time source is, 



H-% + ZH- 1 ^ + 



k 2 gi 
H 2 a 2 



H 



2i / k \3 
"3 fej 



ln(a) +0(c 



(90) 



The solution at late times is straightforward, 



9i(v,k) 



H 



U±-Yln 2 (a) 



W 9 ^HaJ 



0[ ln(a)a 



-3 



(91) 



except for possible homogeneous terms which can be absorbed into a further 
finite field strength renormalization. 



5 Discussion 

We have computed the fully renormalized scalar self-mass-squared at one 
loop (j42M43|) for a general scale factor a(t). For the special case of de Sitter 
(a(t) = e Ht with constant H) our results can be written in a manifestly 
de Sitter invariant form (|46H47|) . Although the computation was simple on 
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account of the conformal invariance of Dirac theory, we have not been able 
to locate a prior result in the literature. 

In any case, our real interest lies in the effect the scalar self- mass-squared 
has on the late time behavior of the scalar mode functions. For that purpose 
we derived the Schwinger-Keldysh effective field equation at one loop 
order. When specialized to de Sitter background we were able show that the 
theory can be renormalized so that there are no significant corrections at late 
times. 

The reason for the null result is that one loop contributions (Fig. 2) 
involve only the conformally invariant part of the theory: the fermion prop- 
agator and the Yukawa coupling. As the introduction explained, significant 
quantum effects during inflation require the participation of at least one 
massless particle which is not conformally invariant. The first instance of 
that for the scalar self-mass-squared would come at two loop order through 
the diagrams of Fig. 3. Evaluating them would be a formidable undertaking, 
but perhaps a possible one in view of the fact that the two loop contribu- 
tions to the scalar self-mass- squared have recently been obtained for massless, 
minimally coupled 4 in de Sitter background [2*3] . 

Recall from the similar analysis of the one loop fermion self-energy 
that the leading one loop correction to the effective field equations came in 
with an extra scale factor compared with the classical terms. It also had 
a factor of f 2 , so we expect the one loop correction to become comparable 
to the classical mode function at roughly f 2 a(t) ~ 1. Two loop corrections 
would be down by an extra factor of f 2 , but should not contain any more 
scale factors. So this is one case in which it is actually reliable to solve the 
effective field equations non-perturbatively with only one loop corrections! 

Of course the physics of the result is that the massless and not conformally 
invariant scalar is catalyzing the production of fermions [13] . The process 
could be throttled if the scalar quickly develops a large enough mass. The 
work done in this paper shows that no such mass occurs at one loop order. 
At two loop order we expect the diagrams of Fig. 3 to induce corrections to 
the effective field equation with an extra factor of a 2 relative to the classical 
terms. However, they will also contain an factor of f 4 , so they would not 
become important until f 4 a 2 ~ 1. That is the same time at which the pro- 
duction of fermions is becoming significant. Because at most one fermion can 
be produced in each state, we expect the end of inflation to see essentially all 
super-horizon fermion modes fully populated. This is a profound difference 
from the usual picture and it seems there must be important consequences. 
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Although we obtained the effective mode equation for an arbitrary scale 
factor a(t), we have only solved it for the special case of de Sitter. For more 
general scale factors there can be stronger asymptotic contributions at late 
times. However, because these vanish for de Sitter, we conclude that any one 
loop late time correction to the mode function must be suppressed by at least 
one factor of a slow roll parameter. The fact that scalar-catalyzed fermion 
production goes to completion ought therefore to be a general feature of slow 
roll inflation. 
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